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Abstract 

Electron-positron annihilation data collected by the BABAR detector near the T(4S) resonance 
are used to study the inclusive decay of B mesons to and D*^ mesons, where the is 
reconstructed using the decay Dg — > cpn . The production fraction of inclusive Ds and the 
corresponding momentum spectra have been determined. The exclusive decays D*~Ds 

are observed with a partial reconstruction technique which uses the soft pion from the D*^ decay 

(*)± 

in association with the reconstructed Dg . The beam energy constraint is used to determine the 
missing mass recoiling against the system, showing a clear signal for this process. From the 
observed rates, preliminary results for the corresponding branching fractions have been obtained. 
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1 Introduction 



The production of the Dg meson in B decays allows one to study the mechanisms leading to the 
creation of a cs quark pair. The main diagram contributing to this decay is shown in Fig. Other 
B decay diagrams also contribute, although at a lower level, but no attempt is made to quantify 
their rate in this paper. As has been pointed out in Ref. the b — ccs decay rate may be large 
and could help to explain the theoretical difficulties in accounting simultaneously for the total 
inclusive B decay rate and the semileptonic branching fraction of the B meson. As a longer term 
goal, the measurement of the rate and momentum spectrum of meson produced in B decays 
beyond the kinematic limit for the process B —i- D*+Xc could be used to study b — > u transitions. 
Despite the fact that purely hadronic final states are more difficult to understand theoretically, 
one may use the particular decay described in this paper to extract Vub/Vcb @- This document 
reports measurements made with the BABAR detector of both the inclusive Dg production rates 
and momentum spectra in B decays and the branching fractions of two specific two-body B decay 
modes involving a Ds meson. The latter measurements are made using a partial reconstruction 
technique. 




d: 



,(*)+ 



Figure 1: The main spectator 
diagram leading to the produc- 
tionofDs mesons in B decays. 



2 The detector and the data sample 

A description of the BABAR detector and the definition of many general analysis procedures can 
be found in an accompanying paper [^. Here only the components of the detector most crucial to 
this analysis are briefly summarized. 

Charged particles are detected and their momenta measured by a combination of a central 
drift chamber (DCH) with a helium-based gas and a five-layer (double-sided) silicon vertex tracker 
(SVT), embedded in a 1.5 T solenoidal field produced by a superconducting magnet. The charged 
particle momentum resolution is approximately {6pt/pt)'^ = (O.OOISpt)^ + (0.005)^, where is in 
GeV/c. The SVT, with typically 10 ^m single-hit resolution, provides vertex information in both 
the transverse plane and in z. 

Particles are identified using a combination of measurements from all the BABAR components. 
Charged particle identification exploits ionization energy loss measured in the DCH and SVT as 
well as Cherenkov radiation measured in a ring imaging detector (DIRC). Electrons and photons 
are identified by the Csl electromagnetic calorimeter. 

Multihadronic events produced in e~^e~ annihilation at the PEP-H collider (SLAC) and collected 
with the BABAR detector have been used in this analysis. These data were taken at the T(4S) 
resonance center of mass energy and at an energy about 40MeV below the BB threshold. The 
integrated luminosity for on resonance data is 7.73 fb~^ and 1.17 fb^^ for off resonance. 
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3 Inclusive production 
3.1 reconstruction 

The mesons are reconstructed using the mode — > (j)7r^ with (p K+K^. In order to 
obtain a sufficiently clean sample, particle identification is necessary. To this end, both energy loss 
(dE/dx) information from the Drift Chamber and the Vertex Detector and the DIRC (a Cherenkov 
imaging detector) are used to identify the kaons produced in the (p decay. 

The selection is based on the likelihoods given by each detector and uses, for each track, the 
ratio of likelihoods for the pion and the kaon mass hypotheses L-j^/Lk- If this ratio is less than 
unity for at least one of the considered subsystems, the particle is selected as a kaon. The DIRC is 
used both in the positive identification mode and the veto mode. A tighter level of identification is 
also available using a total likelihood defined as the product of the likelihoods of each subsystem. 
In this case the track is tagged as a kaon if the ratio of the total likelihoods for the pion and kaon 
mass hypotheses is less than unity. 

Three charged tracks coming from a common vertex are then combined to form a D^ candidate. 
Two oppositely-charged tracks have to be identified as kaons, one of these using the basic criteria 
and the second one using the tighter selection. The K+K~ invariant mass must be within 8MeV/c^ 
of the nominal (j) mass (see Fig. ^). In this particular decay, the cf) meson is polarized longitudinally 
and therefore the angular distribution of the kaons has a cos^ Oh dependence, where the Oh is the 
angle between the K"*" in the (j) rest frame and the cf) direction in the D^ rest frame. We require 
I cos >0.3, thereby keeping 97.5% of the signal while rejecting about 30% of the background. 

Using the selection criteria described above, a reasonably clean signal of D^ is observed (Fig. ^). 
The efficiency averaged over all momenta is (40.5ibl.0)%. It varies as a function of the D^ momen- 
tum and ranges from 30% when the D^ is at rest in the T(4S) rest frame to 55% for p* = 5 GeV/c. 
A clear signal of the Cabibbo-suppressed decay mode D^ — > (pir^ is also observed. A summary 
of the measured signal properties is given in the Table |l|. The final production rates, however, 
are obtained from the invariant mass spectra fitted separately for different momentum intervals 
(Section |3.2D . The measured mass difference mj-ji-m^i agrees with the world average value of 
99.2 ±0.5 MeV/c2 [|. 

Table 1: Fitted parameters for D^ — s- (pi:^ and D*^ — > D^7 decay modes. 





D± ^ (/.7r± 


D:± ^ D3±7 


Fit 


Nds = 18269 ± 202 


N^,± = 3029 ± 151 




M = 1968.5 ±0.1 MeV/c^ 


AM = 143.4 ± 0.3 MeV/c^ 




a = 5.40±0.07MeV/c2 


cjAm = 7.4±0.4MeV/c2 




M^± - M^i = 98.7 ± 0.2MeV/c2 





3.2 Inclusive momentum spectra 

The number of D^ mesons is extracted by fitting the (pii^ invariant mass distribution for different 
momentum ranges in the T(4S) rest frame. The momentum bin width is 200MeV/c, which is much 
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Figure 2: The K+K~ invariant mass spectrum 
for an integrated luminosity of 1.53 fb~^. The 
solid line represents a fit using a Breit-Wigner 
function and a 1st order polynomial. 



Figure 3: The (pw invariant mass spectrum for 
an integrated luminosity of 7.73 fb~^. 



larger than the momentum resolution. The momentum resolution averaged over all momenta 
obtained from the Monte Carlo is 5.6±0.3MeV/c. The fit function is a single Gaussian distribution, 
both for the and the D^. The width of the Gaussians are constrained to be the same and the 
combinatorial background is accounted for by an exponential distribution. The number of in 
the off-resonance data is extracted using the same fit function but with fixed values for M£,± , Mj-,± 
and a obtained from the fit to the on-resonance data. 

The number of reconstructed as a function of their momentum in the T(4S) rest frame is 
shown in Fig. |^ for on- and off-resonance data. The efficiency-corrected momentum spectrum is 
shown in Fig. |5[ 

Table 2: Analytical expressions for the fragmentation functions. 



Name of function Analytical expression 



Peterson et al: f{x,) = ^ (l - J- - ^) 

Collins and Spiller: f{xp) = N + ^e) (1 + (l - ^ 



=^ 

-2 



Kartvelishvili et al: f{xp) = NXp(l 



In order to determine the momentum spectrum from the continuum, on-resonance data 
with momentum higher than 2.45GeV/c and off-resonance data scaled according to the luminosity 
ratio have been fitted after efficiency correction using 3 different fragmentation functions (see 
Table ^). The product of branching fraction, B(D^ — > c/'tt^), times cross-section for production 
from continuum, a{e~^e~ — > D^X), is obtained by integrating the function obtained from the fit 
(Fig. ^). The extracted values and x^/dof bom the fits are shown in Table |^. 
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Figure 4: The momentum spectrum for on- 
resonance data (solid circles) and for scaled off- 
resonance data (open circles) before efficiency cor- 
rection. 



Figure 5: The on-resonance (solid circles) and 
scaled off-resonance (open circles) D5 momentum 
spectrum efficiency-corrected. The solid line is the 
result of the fit using Peterson fragmentation func- 
tion described in text. 



Table ^ shows the contribution of the different sources to the total systematic error for a{e~^e —^ 
D^X) ■ B(Df^ 071"^) • Using the best fit, which is obtained with the Peterson function, we find 
o-(e+e" ^ D±X)-0(D± ^ cpn^) = 8.29±0.41±0.69 pb. From a comparison of the results obtained 
using the other two parameterizations, we assign a conservative systematic error of 2% due to the 
assumed functional form. 

The measured values are in good agreement with previously published results . The momen- 
tum spectrum of the produced in B decays is obtained by subtracting bin-by-bin the value of 
the fit function to the on-resonance data after efficiency correction (Fig. ^). 



Table 3: The parameters for the different fragmentation functions obtained from the fit and the 
measured cross section cr(e"'"e~ — > D^X) • B(D^ — > </)7r^), and the / dof oi the fit. Only the 
statistical errors are given. 



Name of function 


Shape parameter a{e~^e 


BfX) ■ B{Bt - 


^(/>7r±),pb x'^/dof 


Peterson et al: 


e=(12.5±0.6) X 10-2 


8.29±0.41 


1.286 


Collins and Spiller: 


e=(37.6±2.8) x 10^^ 


8.69±0.46 


3.559 


Kartvelishvili et al: 


a=1.91±0.07 


8.63±0.33 


5.338 
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Table 4: The systematic errors for a{e~^e D^X) ■ B{T)^ (/>7r^). 



Source 


Error (%) 


B{(t) ^K+K-) 


1.6 


Particle id efficiency 


0.8 


Tracking efficiency 


7.5 


Luminosity 


3.0 


Total systematic error 


8.3 



3.3 Inclusive branching fraction in B decays 

By integrating the efficiency corrected momentum distribution, a total yield from B meson 
decays of 37050 it 950 events is found. This corresponds to the inclusive branching fraction of 

B{B ^ D±X) = 

where the first error is statistical, the second is systematic and the third is the contribution of the 
(pTT^ branching fraction uncertainty [^. Recognizing that this last uncertainty is common 
to all measurements, our result is slightly higher than the world average (lO.OibO.6 [^) and in good 
agreement with the most precise measurement performed by CLEO The different sources of 
systematic errors are given in detail in Table ^. The dominant uncertainty comes from knowledge 
of the tracking efficiency, which is still the subject of detailed study 

As a cross check of the continuum subtraction procedure, we also subtracted directly the off- 
resonance data scaled by the luminosity ratio for on- and off-resonance. By this means, one obtains 
an inclusive branching fraction B{B T)^X) = 12.0 ± 0.5 ± 1.1%, in agreement with the value 
reported above. 



4 Inclusive D*^ production 

4.1 D*^ reconstruction 

D*^ mesons are reconstructed using the decay D*^ — > D^7 with the subsequent decay — > (pn. 

candidates are selected by requiring the (/>7r invariant mass to be within 2.5 standard deviations 
(it) of the peak value. These are then combined with "single photons" from the event. The later 
are defined by the following criteria: 

• Ey > 50 MeV where is the photon energy in the laboratory frame 

• E* > 110 MeV where E* is the photon energy in the T(4S) frame 



(11.90 ± 0.30 ± 1.07) 



3.6±0.( 



B{Dt 



%, 



(1) 
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Table 6: Systematic errors for B{B D*^X). 

Table 5: Systematic errors for 13{B D±X). 





Error (Vn) 


Source 


H TTTW I /r\ \ 

rji. 1 ui \/0 ) 


Sonrrp 










Signal shape 
Continuum subtraction 


5.0 


Signal shape 


0.9 


1.2 


Background shape 


0.4 


Monte Carlo statistics 


A R 


Continuum subtraction 


1.8 


Bin width 


O.U 


Monte Carlo statistics 


2.0 






Bin width 




Total for yield 


7 7 


0.7 










O.U 


Total for D5 yield 


2.9 




3.6 




2.7 


6(4) K+K-) 


1.6 


Photon efficiency 


2.5 










1.6 


Particle id efficiency 


0.8 








Particle id efficiency 


0.8 


Tracking efficiency 


7.5 




Total systematic error 




Tracking efficiency 
Total systematic error 


7.5 


9.0 


12.0 



• In order to reduce the combinatoric background, the candidate photon should not form a 
■K^ with E*^ > 200 MeV when combined with any other photon in the event. The tt" mass 
window is 115 < My^ < 155 MeV/ c^. 

The distribution of the mass difference AM = Mj-,±^ — Mj^± is shown in the Fig. ^. A clear 
peak with 3030 it 150 events is observed. The parameters obtained from the fit are summarized in 
Table |l[ 

4.2 Inclusive D*^ momentum spectra 

The decay D*± ^ D±7, D± ^ (/>7r^ is used for the measurement of the D*^ inclusive branching 
fraction and the momentum spectrum. The number of D*^ mesons is extracted by fitting the 
AM = Mog-y — Mdj. invariant mass distribution for the different momentum ranges in the T(4S) 
rest frame. A momentum bin width of 400 MeV/ c was chosen. 

The efficiency corrected momentum spectrum is shown in Fig. ^. Both on- and off-resonance 
points corresponding to D*^ mesons produced from the continuum have been fit using different 
fragmentation functions (Table ^). The cross section for D*^ produced from continuum and the 
values of the fit parameters are shown in Table 0. 

Fig. ^ shows the momentum spectrum of D*^ produced in B decays where the Peterson frag- 
mentation function is used for continuum extrapolation. Using this distribution, we find for the 
continuum cross section a{e+e- D*±X) • B{T)f (/>7r±) = 3.48±0.39±0.38 pb. 
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Figure 6: The Dg momentum spectrum efBciency- 
corrected after subtraction of the value of the fit- 
ted curve. The Peterson fragmentation function was 
used for the fit of the continuum. 



Figure 7: AM — Md^^ — A/d^ mass spectrum for an 
integrated luminosity of 7.73 fb~^. The fit function 
is a single Gaussian for the signal and a third-order 
polynomial for the background. 



4.3 Inclusive 0*=*= branching fraction in B decays 

In the same way as for the result, we integrate the efficiency corrected D*^ distribution and 
obtain a total yield from B meson decays of 19300 it 1900 events. From this we find the inclusive 
branching fraction to be 



B{B b;^x) 



3.6 ± 0.9% 
(6.8 ±0.7 ±0.8) X ^ 



%, 



(2) 



where the systematic errors are given in detail in Table 



Table 7: The parameters for the different fragmentation functions, the measured cross section 
iT(e^e" —>■ D*^X) ■ B{J)f- — > </>vr^), and the x^/c^o/ obtained from the fit. Only the statistical errors 
are given. 

Name of function Parameter a{e^e~ —>■ D*^X) ■ B{D^ (t^T^^)-, pb / dof 

Peterson a/.: e=(7.9±0.8)x lO'^ 3.48±0.39 1.260 

Collins and Spiller: e=(19.3±2.3)x 10"^ 3.75±0.42 1.288 

Kartvelishvili et al: q=2.6±0.2 3.61±0.29 1.725 
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Figure 8: The on-resonance (solid circles) and 
scaled off-resonance (open circles) D*^ momentum 
spectrum after efficiency correction. The solid line 
shows the fit using the Peterson fragmentation. 



Figure 9: The D*^ spectrum after efficiency cor- 
rection and continuum subtraction using the result 
of the fit. The Peterson fragmentation function is 
used for the fit of the continuum. 



5 Branching fraction for B° — > D* 0^*^+ decays 

In addition to the measurements of inclusive production rates for and D*^, we have extracted the 
branching ratios for the decays B" — > D*~D+ and B'^ — > D*~D*"'" based on a partial reconstruction 
method. 



5.1 The partial reconstruction method 

As discussed in the introduction, no attempt is made to reconstruct the decays. One combines 
a pion with the reconstructed Ds where the total Ds — vr charge is zero and, assuming that 
their origin is a B^ meson, we calculate the missing invariant mass. This should be the mass 
if the hypothesis is correct^. Without the constraint of the D'' mass, the direction of the B meson 
is unknown. Although its angle with respect to direction can be deduced, the angle (j) around 
this direction is undetermined. Using the beam energy constraint, the missing mass, which still 
depending on the unknown angle (j) of the momentum vector, is computed from: 

mmiss = (£^bcam " Ej^± " K)'^ " {PB " Pd± ~ ' i^) 

In this analysis the missing mass is defined using an arbitrary choice for the angle (p. We use the 
convention that the direction of the B*^ meson lies in the plane {Pn , Pp,{*)±} ■ 



*A11 calculations in this section are performed in the T(4S) rest frame. 
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5.2 Signal extraction 

Fully reconstructed and D*^ are selected by requiring the measured (pir^ mass or Am 



ITT =^7 



-m^^± to be within 2.5 a of the fitted mean value. Because of high combinatorial background 
in the mode with a D*^, one may find several D*^ candidates in an event. Therefore, we form a 
for each candidate defined by: 



X 



Mr 



+ 



2 /M 



rec 
Am 



7\ /fmean 



(4) 



and take the candidate with the lowest value. The Ds —it pairs satisfying the kinematic con- 
straints for the decay B — > D*^Ds are fitted to a common vertex. To reduce further the 
continuum background, we use the event shape variable R2, defined as the ratio of the second to 
zeroth order Fox- Wolfram moment, and require R2 < 0.35. 

The missing mass distributions for the — vr and D*^ — vr are shown in Fig. 10 and 11 
respectively. A clear signal is observed for both decays. The missing mass distribution is fitted 
with the sum of a Gaussian distribution for the signal and a background function given by 

Ci (xo - x\ 

fBix) = S ^ , (5) 



where x is the calculated missing mass, Ci are the parameters of the fit and xq is the end point, 
mD* — mTT = 1.871 GeV/ c^. The results of the fits for both decay modes are summarized in Table |8|. 



> 



200 



150 - 



BaBar 



100 



50 




1.8 



1.825 1.85 



1.875 



M .„ GeV/c^ 

recoil' 



> 




1.85 1.875 
M .„ GeV/c^ 

recoil' 



Figure 10: The missing mass distribution for the 
D^-TT system. The solid line shows the result of the 
fit using the function described in the text. 



Figure 11: The missing mass distribution for the 
D*='=-7r system. The solid line shows the result of 
the fit using the function described in the text. 
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Table 8: Selection criteria and fit parameters for the missing mass distribution in partially recon- 
structed BO ^ D*-Di*^+ decays. 



B D+D*- 


B D*+D*- 


xVdof=l.25 
A^ev = 628±55 events 
mmiss = 1866.7 ± 0.2MeV/c2 
a = 2.31 ±0.15 MeV/c2 


xVdof= 1.09 
A^ev = 195 ± 29 events 
mmiss = 1866.3 ± 0.2MeV/c2 
cr = 2.66±0.36MeV/c2 



5.3 Branching fractions for B° D*-D+ and B° D*-D*+ 

A Monte Carlo simulation of the B*^ D*~Di*^^ decay modes has been used to find the efficiencies. 
It is important to note that the B*^ — > D*~D*^ decay mode contributes to the missing mass 
distribution for the — vr system, even though there is a missing photon from the D*^. We show 
in Table ^ the reconstruction efficiencies for the different modes. 

Table 9: The efficiencies for the partially reconstructed B*^ D*~Ds decay modes. The columns 
show the contribution of the different generated modes to the — vr and D*^ — vr missing mass 
distributions in the signal region. 



True mode 


Reconstructed mode 


D±-7r 




BO ^ D*-D+ 

B° ^ D*-D*+long.pol. 

B° D*-D*+transv.pol. 


32.8±1.8% 
15.8±1.2% 
14.2±1.1% 


9.1 ±0.9% 
6.0 ±0.7% 



Although the amount of feed through from B° D*~D|+ to B° D*~D+ depends on the B*^ 
polarization, one sees from Table ^ that the variation is small. However the reconstruction efficiency 
for B*^ D*~D*"'" has a much larger dependence on the polarization. Since this polarization is 
not known, we use the average efficiency of (7.5 ± 1.5)%, where the systematic error of 1.5% is 
derived by comparing the efficiencies from the two polarization states. This is combined with the 
other systematic errors which are in common with the inclusive branching fractions presented in 
the previous sections. 

For the measurement of the B*^ D*^D*^ branching fraction, the contribution to the missing 
mass peak from B° D*^D+, where a random 7 is associated to the D^, is negligible. The 
contribution of B*^ — > D*~D*^ to the mode B*^ — > D*~D^ is then subtracted to determine the 
branching fraction for the latter. The results are given in Table The first error is statistical, 
the third reflects the uncertainty due to the error in the branching ratio for — > ipir^, and the 
second error represents all remaining systematics. This last is dominated by the uncertainty due 
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to the dependence of the efficiency on the polarization of the final state. 



Table 10: The measured branching fraction for B° D* D+ and B° D* D*+. 

Preliminary 



B D+D*- 


B D*+D*- 


B= (7.1±2.4±2.5±1.8) x lO^^ 
PDG: S=(9.6±3.4) x lO^^ 


e=2.54±0.38±0.53±0.64% 
PDG: B=2.0±0.7% 



Finally, one should note that the reconstructed B D*^D*+ events should allow us to measure 
the polarization of the D*^ in these decays and therefore, in future analyses, it will be possible to 
reduce the systematic error from this source. 



5.4 Background cross checks 

In order to investigate further the shape of the background which is subtracted for estimating the 
signal, we have compared the Monte Carlo to the data. Several types of backgrounds contribute in 
the signal region: 

1. Fake D^*^^ and a random pion (for example coming from the other B). 

2. Fake Ds*^^ and correlated pion (for example coming from the same B). 

3. True Ds and a random pion. 

4. True D^*'*^ and a correlated pion. 

Table O shows the different types of backgrounds and the methods which are used to determine 
their level. Background types 1 + 3 are obtained by flipping the Dg direction. Background 
types 1 + 2 are extracted using the sidebands of the Ds mass distribution. For this purpose, 
we take 1.89< M^i <1.95 and 1.985< Mq± <2.05 GeV/c^ for the D^-vr system, and AM^*± 

170< AMj-,.± <300 MeV/c^ for D*=^-7r. By ffipping the D^*^"^ direction for the sidebands we find 
the contribution of background type 1. Therefore the difference between the distributions for ffipped 
and non-flipped Dg direction for the sidebands gives the type 2 background contribution and 



thus it is possible to find the contribution of background types 1 + 2 + 3 from data alone. Fig. 12 
and |l^ show the resulting signal after their subtraction. The remaining background component is 
quite small and is estimated from the Monte Carlo. To ensure that the simulation reproduce the 
data well, a systematic comparison is made for the missing mass distribution obtained from the 
D^ signal region, the D^ sideband region, and the wrong-sign D^ — vr combinations both in the 
D^ signal and the D^ sideband regions. The ratio (Data- Monte Carlo)/Monte Carlo for all these 
cases are determined as a function of the missing mass. We find good agreement within the errors 
in all cases. Table ^ summarizes this result by showing the ratio integrated over the missing mass 
region 1.78 to 1.87GeV/c^ for all distributions except that with the signal, for which the range 1.78 
to 1.85GeV/c2 is used. 
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Table 11: The different data samples which can be used to determine the background in the D'^ 
signal region. 



Background 


Flip D^*^^ D^*^^ Side-bands Side-bands flip D^*^^ 


1. Fake Dg + random vr 

2. Fake Dg + correlated vr 

3. True Dg + random tt 

4. True Dg + correlated tt 


XX X 
X 

X 



Table 12: The comparison of the different data samples with Monte Carlo. 







B° D*-D*+ 


Sample type 


(A^data - iVMc)/A^MC 


xVdof 


(A^data - iVMc)/A^MC 


xVdof 


D± Signal 


0.051±0.025 


1.008 


0.103±0.057 


1.058 


Flip D± 


-0.043±0.031 


0.841 


-0.041±0.064 


0.832 


D± Sideband 


0.006±0.018 


1.391 


-0.031±0.053 


1.194 


Flip Sideband 


0.015±0.021 


1.627 


0.084±0.069 


1.690 


Wrong Sign 


-0.031±0.029 


0.987 


0.010±0.063 


1.088 


WS, D± Sideband 


0.030±0.020 


1.311 


0.034±0.065 


1.487 



6 Conclusion 



The production of D^*^* at the T(4S) energy (and 40 MeV below) has been studied with the BABAR 
detector. Preliminary measurements of branching fractions for t'inclusive production and for the 
exclusive decays B° D*~D*"'" have been performed. The following cross sections have been found 
for production in the continuum: 

a{e+e- Bfx) ■ B{Bf c^tt^) = 8.29 ± 0.41 ± 0.69 pb , 

(7(e+e~ ^ D*^X) • B(Dt </'7r^) = 3.48 ± 0.39 ± 0.38 pb . 
Using the on-resonance data, the inclusive branching fraction for the B meson decays 

3.6 ±0.9^ 



B{B ^ D±X) = 



B{B D*±X) 



(11.90 ±0.30 ±1.07) X 
(6.8 ±0.7 ±0.8) X 



B(Dt - 

3.G ± 0.9% 



TT 



% 
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Figure 12: The missing mass distribution for tlie 
D^-TT system from data (points) and Monte Carlo 
(histogram) after background subtraction (see text). 



Figure 13: The missing mass distribution for 0*"*=- 
TT system from data (points) and Monte Carlo (his- 
togram) after background subtraction (see text). 



have been measured. Finally the decays B*' D* D+ and B'^ — > D* D^*'*^ have been observed 
using a partial reconstruction technique and the following branching fractions have been determined: 



B{B- 



D+D* 



') = (7.1±2.4±2.5±1.8) X 10"^ ; 
D*+D*-) = (2.54 ± 0.38 ± 0.53 ± 0.64)% . 



The results obtained are in a good agreement with previous measurements by other experiments. 
The measurement of inclusive branching fraction of D*^ from B decay has been obtained for the 
first time. 
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